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We study one dimensional quasiperiodic system obtained from the tight-binding model on the 
square lattice in a uniform magnetic field with the spin orbit interaction. The phase diagram with 
respect to the Harper coupling and the Rashba coupling are proposed from a number of numerical 
studies including a multifractal analysis. There are four phases, I, II, III, and IV in this order from 
weak to strong Harper coupling. In the weak coupling phase I all the wave functions are extended, in 
the intermediate coupling phases II and III mobility edges exist, and accordingly both localized and 
extended wave functions exist, and in the strong Harper coupling phase IV all the wave functions 
are localized. Phase I and Phase IV are related by the duality, and phases II and III are related 
by the duality, as well. A localized wave function is related to an extended wave function by the 
duality, and vice versa. The boundary between phases II and III is the self-dual line on which all 
the wave functions are critical. 

In the present model the duality does not lead to pure spectra in contrast to the case of Harper 
equation. 

PACS numbers: 71.23.An, 71.23.Ft 



I. INTRODUCTION 

The quasiperiodic systems are intermediate between 
periodic and random. Electronic wave functions are al- 
ways extended in periodic systems. In one dimension, 
wave functions are always localized in random systems. 
The one-dimentional quasiperiodic systems show various 
interesting wave functions and energy spectral characters 
(for a review, see [![). 

Consider the one-dimensional tight-binding model: 

- 1p n+ l - 4> n -i + V(um)lp n = Elp n , (1) 

where V(x) is a periodic function with period 1, i. e. 
V(x + 1) = V(x). When lj is an irrational number, this 
becomes quaisperiodic. In particular, the model defined 
by 

V{x) = 2Acos(2ttx), (2) 

is known as Harper equation. It has Aubry duality [2j 
which relates the region with A > 1 and the region with 
A < 1. Thus A = 1 is the self-dual point. From this dual- 
ity and Thouless's formula for Lyapunov exponent [3j, it 
has been argued that all the wave functions are extended 
for A < 1, while all the wave functions are localized for 
A > 1. At A = 1, all the wave functions are critical 
Thus extended, localized, and critical wave functions do 
not coexist in the energy spectrum as it is pure at each 
value of A; the spectrum for A < 1 is purely absolutely 
continuous, it is purely singular continuous for A = 1, 
and it is purely dense point for A > 1. 

A pure spectrum, however, is not a general feature of 
the one-dimensional quasiperiodic systems. It is known 
that the model in Ref. 4 which has 

V(x) — Atanh(/xcos (2irx))/ tanh/i, (3) 



Q and the genenalized Harper model Q in which V(x) 
has two periods: 

V(x) = Ai cos(27ra;) + A m cos(2m7ra;), (4) 

have non-pure spectra and mobility edge structures. 
If we choose 

f-A, m -a G <x<m 

I +A, m < x < m + 1 — <tq, 

where oj ~ ao = (v5 — l)/2 (the inverse of the golden 
mean) and m is an arbitraly integer, the potential follows 
the Fibonacci sequence which is constructed recursively 
as 

{M e+1 } = {M e _,}{M ( }, (6) 

with {M } = B, {Ah} = A. Thus {M 2 } = BA, {M 3 } = 
ABA, {Mj} = BAABA, ■■■ . This is called the Fibonacci 
model[9|, [lOj] . From the renormalization group study it is 
shown that the energy spectrum is purely singular con- 
tinuous and all the wave functions are criticalpj], H3 |. 
This property is proven rigorously [13J- 

The spin Hall effect is induced by spin orbit coupling 
in which spin current is in the transverse direction to an 
applied electric field in the absence of a magnetic field 
[lj]. This phenomenon has been attracting much inter- 
est because of rich physics and potential industrial ap- 
plications. This is in contrast to the quantum Hall effect 
in which an applied magnetic field breaks time reversal 
symmetry. But the theory of quantum Hall effect plays 
an important role in some aspects [l5|, [l6|, [l?], [HI • 

Some of the basic models of the spin Hall effect involve 
the Rashba spin orbit coupling. It can change properties 
of electron systems essentially. For example, the level- 
avoiding gaps for single electron system with disorder on 
a cylinder shows a new kind of universal behavior in the 
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presence of the Rashba spin orbit coupling. The distribu- 
tion shows an exponential decay rather than the Gaussian 
one, which can not be explained by the Gaussian unitary 
ensemble [l9| . 



In this paper the tight binding model on the square lat- 
tice in a magnetic field is considered. The Rashba spin 
orbit interaction is also included. The model is reduced to 
Id and it is quasiperiodic when magnetic flux per plaque- 
tte is irrational. It is a generalization of Harper equation 
and we study the properties of the energy spectrum and 
the wave functions of this Id quasiperiodic system. We 
find four phases: the weak coupling phase I, the strong 
coupling phase IV, and the intermediate coupling phases 
II, III (see Fig. [T]). All the wave functions are extended 
in the weak coupling phase I, and all the wave functions 
are localized in the strong coupling phase IV. In contrast 
to Harper equation localized and extended wave func- 
tions coexist, and accordingly mobility edges appear in 
the intermediate coupling phases II and III. The bound- 
ary of these phases is the self-dual line where all the wave 
functions are critical, and the spectrum is singular con- 
tinuous. The self-dual line is a critical point of a phase 
transition between a localized wave function and a ex- 
tended wave function. 



The Rashba spin orbit interaction does not break 
Aubry duality, but introduces a coupling between two 
Harper equations (one for spin up and the other for spin 
down). Our result suggests that an extended state in 
phase I could be either localized or extended in phases 
II and III. This is also true for a localized state in phase 
IV. Thus the spin orbit coupling enhances tendency to- 
ward neither localization nor derealization. This sug- 
gests that the existence of mixed spectra is a rather com- 
mon property of the Id quasiperiodic system. The dual- 
ity itself does not guarantee a pure spectrum. 



The organization of this paper is as follows: in Sec. |TT] 
we introduce two dimensional elecrons with the Rashba 
spin orbit interaction on the square lattice in a uniform 
magnetic field and obtain the one-dimensional quasiperi- 
odic system. The duality which is a generalization of 
Aubry duality of Harper equation is discussed. The pro- 
posed phase diagram is shown and explained in Sec. IIIII 
In Sec. HVI the scaling approach is explained. The analy- 
sis of the total band widths is presented which supports 
the proposed phase diagram. In Sec. [V] the multifrac- 
tal method for an energy spectrum and wave functions 
is introduced. The scheme for identifying phases with 
absolutely continuous, singular continuous, dense point, 
and mixed energy spectra are given. In Sec. IVII exam- 
ples of extended, localized, and critical wave functions 
are shown. 



II. HAMILTONIAN AND DUALITY 

A. Electron hopping on the square lattice in a 
magnetic field with the spin orbit interaction 

Consider tight binding electrons on the square lattice 
in a perpendicular magnetic field with the Rashba spin 
orbit interaction [l4j]. The problem involves a Hamilto- 
nian: 

H - -Vft c* c e^+ic 1 c e l0V ^A 

11 ~ / A xL "n+l,m L "i'i,mZ T L y'-n,m+l'-ri,m^ ) 

+ H.c, (7) 

where c„, m = 4 (c», m T, Cn,mi) and c\ m = {c\ l mV c ] n mi ) 
(t is a transpose). The transfer integrals in x- and y- 
directions are given by t x and t y , respectively, and Xr is 
the Rashba spin orbit coupling. When = 0, Eq. ([7]) 
reduces to the Rashba model on the square lattice. 
We take the Landau gauge: 

0n, m = 0, %, m = 2 ^<K (8) 

and flux per plaquette is given by 2ir(j). With this choice 
of gauge the Hamiltonian has translational symmetry in 
y-direction and one can write 

W)= E e^n^JO), (9) 

n,m,a=(1 ,1) 

then H\ip) = E\ip) leads to 

+ A fl {ipn+1,1 - ipn-i,l - 2sin(27r0n - k v )ip n>if } = E<ip n j, 

-t x (ip n+1< i +-0„-u) -2t v cos {2it(jm - k y )ipn t i 

+ X R {-ip n+ ij + ipn-i,l ~ 2sin(27n/>n - k v )^ n ^} = Eip nt ±, 

(io) ■ 

Let us make a gauge transformation by exchanging x- 
and y-directions. The new gauge is 

9f hm = 2^m, flv m = 0, (11) 

then we have 

- £y(/m + l,T + /m-l,t) - COs(27T0TO - k x )f m ^ 

+ Afl {/m+i,i - fm-i,l ~ 2sin(27r^m - k x )f md } = Ef m ^, 

- ty(f m +l,l + /m-i,;) - 2t x C0s(27T(/>m - k x )fm,l 

+ Afl {-/m+i.T + /m-i.T - 2sin(27T0m - k x )f mA } = Ef m ^, 

(12) 

where, 

W= J- e^„ h „clJO), (13) 

n,m, <t=(T,1) 
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The Id wave functions in Eqs. (jTOj) and IT2|) are related 
by a Fourier transformation, 

Vn,a = a(a)e- ik * n ^ e l(2 ^"-^ )m / m , ff , (14) 
m 

where 

a(T) = l, a(l)=i. (15) 

B. Quasiperiodic system in one dimension and the 
duality 

Equation (|10[) represents a tight binding model in Id. 
When is a rational number, <f> = p/q, where p and q 
are mutually prime integers, it is periodic with period 
q. Thus Bloch's theorem is valid and the wave functions 
are of Bloch form and extended. The energy spectrum 
consists of 2q bands. 

When 4> is an irrational number, which could be con- 
sidered as p — > oo and q — * oo, Bloch's theorem cannot 
be applied. Thus it is not certain whether the wave func- 
tions are extended or not. In addition, since the number 
of bands 2q is infinite, the energy spectrum is expected 
to be a Cantor set. This is the localization problem we 
study. 

In order to construct a phase diagram the duality plays 
a key role. Compare Eqs. I|10p and (fT!?|) and one notices 
that they have the same form except that t x and t y , as 
well as k x and k v , are exchanged. The wave functions 
are related by the Fourier transformation (|14p . Let us 
introduce a parameter, 

\h = % (16) 

then the system is self-dual if Xh = 1. A wave function 
for Xh = a is related to the one for Xh = 1/a by the 
Fourier transformation (|14p . Thus, if a wave function 
is exponentially localized for Xh = a, for example, the 
corresponding wave function for Xh = 1/a is extended. 

If Xr = 0, there is no coupling between electrons with 
opposite spins and the system is reduced to Harper equa- 
tion. In this case the duality is called Aubry duality [2j. 
Thus the duality in the present system is a generalization 
of Aubry duality. 

III. PHASE DIAGRAM 

If a energy spectrum is pure, there are three possi- 
bilities: absolutely continuous, dense point, and singu- 
lar continuous. For an absolutely continuous spectrum 
all the wave functions are extended, for a dense point 
spectrum all the wave function are localized, and for a 
singular continuous spectrum all the wave functions are 
critical. Another possibility is a mixed spectrum which 
consists of absolutely continuous parts and dense point 




FIG. 1: (Color online) Proposed phase diagram. Phase I 
and Phase IV are dual each other, as well as Phase II and 
Phase III. The system is self-dual on the red line, Ah = 1; 
the energy spectrum is purely singular continuous and all the 
wave functions are critical. Phase I: The energy spectrum 
purely absolutely continuous and all the wave functions are 
extended. Phase II and Phase III: The energy spectrum is 
mixed and consists of absolutely continuous parts and dense 
point parts. Phase IV: The energy spectrum is purely dense 
point and all the wave functions are localized. 

parts. Therefore extended wave functions and localized 
wave functions coexist in an energy spectrum in this case. 

In Fig. [1] we show the proposed phase diagram for the 
present quasiperiodic Id system. In the week coupling 
phase I the spectrum is absolutely continuous and all the 
wave functions are extended. In the intermediate phases 
II and III the spectrum is mixed; the absolutely contin- 
uous parts and dense point parts coexist and they are 
separated by mobility edges. These two phases are re- 
lated by the duality. So the absolutely continuous parts 
and dense point parts of the spectrum are exchanged by 
the duality which is discussed in Sec. |TT]B. The extended 
wave functions and the localized wave functions are re- 
lated by the transformation (fT4]) . 

The boundary of Phase II and Phase III is the self-dual 
line, Xh = 1- We expect that all the wave functions are 
critical and the spectrum is purely singular continuous. 

In the strong coupling phase IV the spectrum is dense 
point and all the wave functions are localized. This 
strong coupling phase IV and the weak coupling phase 
I are related by the duality. The corresponding wave 
functions are related by the transformation (fT4"|) and all 
the wave functions are extended in Phase I. 



IV. THE SCALING APPROACH 

We analyze the energy spectrum and the wave func- 
tions by the scaling approach. We choose the inverse of 
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the golden mean, oq ■ 
4>. Note the relation, 



(y/E— l)/2, for irrational number 



lim 



F, 



Ft 



(17) 



where the Fibonacci numbers Fi are defined recursively 
by F e+ i = Fe-i + F e , with Fo = F 1 = 1. We take ra- 
tional approximations <pi — Fe-i/Fg in which the system 
is periodic with period Ft. Then systematic studies for 
large I are performed. 

When <j> £ = Fz-i/Ft = p/q, Eq. (JTOj) is 
periodic with period q. From Floquet's or 

Bloch's theorem we have ip n + q = e, lkq ^Pn {n = 
1,2, •••,<?)■ Thus, if we introduce a vector 

"0 = (V>l,T>V'U>V'2,'bV'2 ) .l,--- j^-l.T'^-ia^g.T'V^-l), 

Eq. (flTJ|) is reduced to an eigenvalue problem of a 2g x 2g 
matrix: 
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(18) 



where 



-An — 



, } , t y cos{2 / K(j>in) Xn sin (2ir4>en) 
k \r sin (2ircf>gn) t y cos (27n^n) 



B = 



-Ah 



and 



C = 



Xr 



Xr 



(19) 



(20) 



(21) 



We set fc y = for simplicity in Eq. (TiTJ)) . The first 
Brillouin zone of this Id system is 

--<**<-. (22) 

Then, the energy spectrum has 2q = 2Fg bands. 

A. Total band width and the phase diagram 

Quantities which are useful in investigating phase tran- 
sitions of quasiperiodic systems are sums of band widths 
of rational approximants which are denoted by i30, 0j- 
If B is non-zero in the quasiperiodic limit, I — ► oo, the 
energy spectrum has absolutely continuous parts and the 
corresponding wave functions are extended. On the other 
hand, wave functions are localilzed or critical when B 
tends to zero. 
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FIG. 2: Total band widths of Harper equation. The limit 
q — > oo represents the quasiperiodic system. 



When Xn = 0, there is no spin orbit interaction and 
the problem is for electrons on the square lattice in a 
magnetic field, namely Harper equation. As a canonical 
example of the scaling approach, results for B in Harper 
system are shown in Fig. [2j For A# < 1, B approach 
a finite value for q — > oo while for Xh > 1, B tends 
to zero rapidly. So the energy spectrum is dense point 
for Xh > I- The result for A# < 1 suggests that the 
energy spectrum contains an absolutely continuous parts. 
But with the aide of the duality which connects the two 
regions it can be concluded that the energy spectrum is 
purely absolutely continuous for Xh < 1 ■ For Xh = 1 the 
total band width behaves as 



B~q- 



(23) 



This behavior suggests a singular continuous spectrum 
and the exponent S is an index which characterizes the 
criticality and is estimated in Ref. 4 to be 1.00. 

The scaling method for B is applied to the present 
model Eq. (TIT)]) . In Fig.[TJ the boundary between Phase 
III and Phase IV is denoted by g(Xn)- Examples of data 
which determine g(Xn) is shown in Fig. [3Ja). For > 
g(Xn) B tends to a finite value, while for A# < g{Xjj) B 
tends to zero. In Pase III there are absolutely continuous 
parts in the energy spectrum, but there is no absolutely 
continuous parts in Phase IV. 

The green error bars in Fig. [T] are determined by ob- 
serving that, just outside the errors bars, the total band 
widths tend to a finite value or go to zero exponentially 
within the numerical accuracy. 

The boundary between Phase I and Phase II, g'(A//), 
is related to g(Xn) by the duality as 



1 



H 



(24) 



and is shown in Fig.[TJ In Phase I and Phase II B always 
tends to a finite value as shown in Fig. [3fb). 
The function is estimated as 



g(X H ) = a(X H -if + b(X H 
+c{X H - if- 2 , 



1 



(25) 
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FIG. 3: Total band widths for X R > 0. The limit q -> oo 
represents the quasiperiodic system. 



where a = 0.14 ± 0.02, 6 = -0.72 ± 0.09, c = 1.34 ± 0.09 
and B = 3.03 ± 0.08. 

On line Ah = 1 the system is self-dual. The scaling 
analysis shows that the system is always critical. Exam- 
ples of the scaling analysis are given in Fig. [3^c). The 
scaling index defined by Eq. (j23|) is independent of A.r 
and is given by S — 1.00. 



V. MULTIFRACTAL ANALYSIS 

We further investigate scaling property of the energy 
spectrum and the wave functions by the multifractal 
analysis [2(| to strengthen the proposed phase diagram. 



In Ref. [2tJ a general formulation of the multifractal the- 
ory is introduced which reformulates the theory along the 
line of the standard statistical mechanics. This method 
can be applied to both energy spectra and wave functions 
of quaisiperiodic systems but the formalisms are different 
for the two cases. 

The exact f(a) spectrum is obtained for the six-cycle 
wave function of the Fibonacci model [13, k3| . 



A. Multifractal anaysis of an energy spectrum 

In order to characterize scaling of an energy spectrum 
we define appropriate scaling indices and the entropy 
function. At £th level of rational approximations, we 
have N — 2Fg bands and their widths are denoted by Aj 
(£ = !,••• , N). Let us define a scaling index for Aj by 



1 



■ In A ,- 



We also define an entropy functilon by 
5(e) = lim iln^(e), 



(26) 



(27) 



where Qi{e)de is the number of bands whose scaling index 
lies between e and e + de. As in the formalism of statis- 
tical mechanics, we introduce a "partition function" and 
a "free energy" which are defined by 



N 



and 



F{8) = lim -AnZ t {B). 



(28) 



(29) 



Once the free energy is calculated, the enetropy function 
is obtained by the Legendre transformation: 



S(e) = F({3) + Be, 



dF{B) 
dB ' 



(30) 



Thus, by changing B, one can pick a value of e, and the 
corresponding 5(e) is calculated. 

The index e which represnets scalilng of the Lebesgue 
measure of the energy spectrum can actually be related 
to singular behavior of the density of states. At £th level 
of approximation, each band carries the same number of 
states, 



Pi 



1 



1 

—t 



(31) 



where r is the golden mean: (V5 + l)/2. The total num- 
ber of states is normalized to unity, i.e., YliPi = 1- An 
index on which represents singular behavior of the den- 
sity of states is defined by 



Pi 



Af\ 



(32) 
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FIG. 4: (Color online) f(a) for energy spectrum on the self- 
dual line. At Xr — it corresponds to that of Harper equa- 
tion. 



Since pi = 1/N = r and ~ e , we have a relation: 

at = Iiit. (33) 

The spectrum of singularity introduced by Halsey et al. 
[2fj| is given by 



n'(e,a) = <A)-A») > 



(34) 



where fl'(e, a)deda is the number of bands whose scaling 
index lies between e and e+de, and a and a+da, and (A) 
is a representative value of A, (A) = exp (— £e). Then 
/(cv) is related to the entropy function S(e) by 



/(«) 



S(e) 



(35) 



The point here is that wether the limit £ — > oo in Eq. ([77]) 
exists or not. Only when it exists the energy spectrum is 
multifractal. 

We apply the method above to the energy spectrum on 
the self-dual line, Ajj = 1. It is possible to obtain f(a) 
numerically by extrapolation of large £, and some of the 
results are shown in Fig. [5] These /(cv)'s exhibit smooth 
curves in intervals [a m i n , a max ] . This is an evidence that 
the energy spectrum is singular continuous and we con- 
clude that all the wave functions on the self-dual line is 
critical. 



B. Multifractal analysis of wave functions 

Take a period of the lattice which has a Fibonacci num- 
ber of lattice points. Normalize the length of the period 
by unity by putting the lattice spacing as 



Ft T l 



(36) 



The scaling index for U, which is given by U — exp(— £ti), 
is a constant: 



1 



• In I ~ In r. 



(37) 



We define a probability measure from a wave function by 

Pi = l^l 2 , (38) 
which is normalized as 



F t Ft 

$> = £M a = i. 

i=l i=l 



(39) 



The scaling index for pi is defined by 

Pi = I?, (40) 
and the entropy function for a is defined by 



S(a) = lim — \nQe(a), 

£— >OG t 



(41) 



where Q,f>{a)da is the number of sites which have an index 
between a and a + da. 

As in the previous case of the energy spectrum we in- 
troduce a partition function: 



z l { q ) = Y J pi 



and a free energy: 



G(q)= lim -lnZ e (q). 



(42) 



(43) 



The entropy function is given by the Legendre transfor- 
mation: 



S(a) = G(q) + qae, 

1 dG(q) 
a = . 

e dq 

Function f(a) in this case is defined by 

Q( a ) = r /(a) , 
and is related to the entropy functiton as 
tt \ S(a) 



(44) 



(45) 



(46) 



We calculate f(a) for finite Fibonacci indices £ and 
extrapolate them to £ — > oo. Actually only a part of 
/(a) is required to distinguish localized, extended, and 
critical wave functions. For extended wave functions, the 
probability measure scales as pi ~ 1/N = I, thus we have 
Q^min = 1- For localized wave functions, on the other 
hand, pi is finite (a = 0) only on a finite number of sites 
and is exponentially small (a = oo) on the other sites, 
thus we have f(a) which is defined at a = [/(0) = 0] 
and a = oo [/(oo) = 1]. For critical wave functions, a has 
a distribution, i.e., f(a) is a smooth function defined on a 
finite interval [a m ; n , a m ax][l|- Thus we need to calculate 
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amin which corresponds to the largest |V'j| 2 to distinguish 
extended, localized and critical wave functions. Namely, 

OL m in = 1 for an extended wave function, 

O-min 7^ 0, 1 for a critical wave function, 

OL m in = for a localized wave function. (47) 

At £th approximation of the Cantor set sectrum, there 
are N = 2F? energy bands. We obtain a wave function 
with (k x ,k y ) = (0,0) at each band. Examples of deter- 
mining a m i n for wave functions in Phase I and Phase IV 
are shown in Fig. [5] (a). The wave functions of the low- 
est band are denoted as i = 1. For this wave function 
a min extrapolates to if (Xh,Xr) — (1-05,0.01), and to 
1 if (A H , A fl ) = (1/1.05,0.01) (see Fig. 0(a)). In this 
manner it is determined that a m i n = for all the wave 
functions if (Xh,Xr) = (1.05,0.01), and a m i n = 1 for all 
the wave functions if {Xh,Xr) — (1/1.05,0.01) (see Fig. 
[5](b)). The wave functions at the center of the spectrum 
(i = Fg) give the same results for a m i n . In this way the 
multifractal analysis of wave function also shows that all 
the wave functions are extended in Phase I and all the 
wave functions are localized in Phase IV. 

Examples of the results for Phase II and Phase III 
are shown in Fig. [SJ For the wave function of the low- 
est band (i — 1) a min extrapolates to if (A#,Ar) = 
(1.05,0.8), and to 1 if {X h ,Xr) = (1/1.05,0.8) (see Fig. 
[6] (a)). On the other hand, a m i n for the wave func- 
tions at the center of the spectrum (i — Fg) a m in ex- 
trapolates to 1 if (\h,Xr) = (1.05,0.8), and to if 
(Ah, Afl) = (1/1.05,0.8) (see Fig. ©(a)). This suggests 
that there exist a number of mobility edges in the energy 
spectrum. The existence of mobility edges is numerically 
shown using this multifractal analysis of wave function, 
and an example of the data are plotted in Fig. [S] (b). 

The duality discussed in Sec. [TT] E> dictates that, if a 
wave function with (Xh,Xr) = {a, (3) is extended, the 
corresponding wave function with (A#, Xr) — (1/a, 0) is 
localized, and vise versa. This behavior is clearly seen in 
Fig. El (b). 
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(b) Xr = 0.01 




o 0.5 1 



i/2F 18 



FIG. 5: (Color online) (a) Plots of a m in vs. l/£ for the wave 
functions of the lowest band (i = 1) and the center band 
(i = Ft), (b) Plots of a m in for the wave function of each 
band. The number of wave functions is 2Fi$ (Fig = 4181). 

of wave functions in Phase IV in which all the wave func- 
tions are localized. Examples of critical wave functions 
on the self-dual line are shown in Fig. [TTJ 



VI. WAVE FUNCTIONS 

In this section examples of the wave functions in pa- 
hases I, II, III, and IV; and on the self-dual line are dis- 
played. 

Example of wave functions in Phase I are shown in 
Fig. [7] All the wave functions are extended in Phase 
I. Examples of wave functions in Phase II are shown in 
Fig. [SI We see that the wave function at the edge of the 
energy spectrum is extended and the wave function at 
the center of the energy spectrum is localized. Examples 
of wave functions in Phase III are shown in Fig. [9] We 
see that the wave function at the edge of the energy spec- 
trum is localized and the wave function at the center of 
the energy spectrum is extended. Figure. ITTJ1 is examples 
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FIG. 6: (Color online) (a) Plot of a m i n vs. l/£ for the wave 
functions of the lowest band (i = 1) and the center band 
(i = i 7 ^). (b) Plots of a m in for the wave function of each 
band. The number of wave functions is 2Fis (Fis = 4181). 
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FIG. 7: Square moduli of wave functions in Phase I. 
(Ajy, Xr) = (1/1.05, 0.01), £= 19, Ft = 6765. (a) At the edge 
of the energy spectrum, (b) at the center of the energy spec- 
trum. 



an 
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2000 



FIG. 8: Square moduli of wave functions in Phase II. 
(Ah, \r) = (1/1.05, 0.8), I = 19, Ft = 6765 (a) At the edge of 
the energy spectrum, (b) at the center of the energy spectrum. 



9 




FIG. 9: Square moduli of wave functions in Phase III. FIG. 11: Square moduli of wave functions on the self-dula 

{Xh,Xr) = (1.05, 0.8), £ = 19, F t = 6765. (a) At the edge line, (Ah, Ah) = (1,0.2), £= 19, F e = 6765. (a) At the edge of 

of the energy spectrum, (b) at the center of the energy spec- the energy spectrum, (b) at the center of the energy spectrum, 
trum. 



(a) 

0.02 




FIG. 10: Square moduli of wave functions in Phase IV. 
(Ah, Ah) = (1.05, 0.01), I = 19, Fig = 6765. (a) At the edge of 
the energy spectrum, (b) at the center of the energy spectrum. 
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